Supersymmetric Optical Structures 
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We show that supersymmetry can provide a versatile platform in synthesizing a new class of op- 
tical structures with desired properties and functionalities. By exploiting the intimate relationship 
between superpatners, one can systematically construct index potentials capable of exhibiting the 
same scattering and guided wave characteristics. In particular, in the Helmholtz regime, we demon- 
strate that one-dimensional supersymmetric pairs display identical reflectivities and transmittivities 
for any angle of incidence. Optical SUSY is then extended to two-dimensional systems where a link 
between specific azimuthal mode subsets is established. Finally we explore supersymmetric photonic 
lattices where discreteness can be utilized to design lossless integrated mode filtering arrangements. 
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Supersymmetry (SUSY) emerged within quantum field 
theory as means to relate fermions and bosons [T]-[6]. 
In this mathematical framework, these seemingly very 
different entities constitute superpartners and can be 
treated on equal footing. Transitions between their re- 
spective states require transformations between commut- 
ing and anti-commuting coordinates-better known as 
supersymmetries. The development of SUSY was also 
meant to resolve questions left unanswered by the stan- 
dard model [71 , such as the origin of mass scales or the 
nature of vacuum energy, and to ultimately link quan- 
tum field theory with cosmology towards a Grand Unified 
Theory. Moreover, SUSY has found numerous applica- 
tions in random matrix theory and disordered systems 
[12 . Even though the experimental validation of SUSY 
is still an ongoing issue, some of its fundamental con- 
cepts have been successfully adapted to non-relativistic 
quantum mechanics (QM). Interestingly, in this context, 
SUSY has led to new methods in relating Hamiltonians 
with similar spectra. In this regard, it has been used to 
identify new families of analytically solvable potentials 
and to enable powerful approximation schemes [8-11 . 
Recently, SUSY schemes have been theoretically explored 
in quantum cascade lasers [13 and ion-trap arrangements 
[14 . Clearly of interest will be to identify other physical 
settings where the rich structure of SUSY can be directly 
observed and fruitfully utilized. 

In quantum mechanics, SUSY establishes a relation- 
ship between superpartners through the factorization of 
an operator, i.e., C^^^ = AM, where f denotes the Her- 
mit ian adjoint. In this respect, the superpartner is de- 
fined through C^'^^ = AA^ ^ from where one finds that 
^£(1) = AA^A = CP)A and A^ [P^ = A^ AA^ = 
C'^^^A^. It then follows that the two eigenvalue prob- 
lems £(1^2)^^(1,2) ^ ^(1,2)^^(1,2) y-^y identical spectra 
^(1) ^ ^(2)^ Moreover, the SUSY operators A^ and 
A pairwise transform the eigenfunctions of the respec- 
tive potentials into one another: X^^^ ex A^X^'^'^ and 
X^^^ (X AX^"*^^ [8 . In addition, supersymmetry demands 
that A annihilates the ground state of £^^\ Therefore 



the corresponding eigenvalue is removed from the spec- 
trum of vC^^^ . If however A does not annihilate the ground 
state of C^^^ , then the two operators share the exact same 
spectrum (including the fundamental state), and SUSY 
is said to be broken. In the language of superpotentials, 
this may also be characterized through the Witten pa- 
rameter [6l E] . 

In this Letter we show that optics can provide a fer- 
tile ground where the ramifications of SUSY can be ex- 
plored and utilized to realize a new class of functional 
structures with desired characteristics. In particular we 
demonstrate that supersymmetry can establish perfect 
phase matching conditions between a great number of 
modes-an outstanding problem in optics. In this vein, we 
illustrate the intriguing possibility for preferential mode- 
filtering where the fundamental mode of a structure can 
be selectively extracted. Moreover, in the Helmholtz 
regime, SUSY endows two very different scatterers with 
identical reflectivities and transmittivities irrespective of 
the angle of incidence. Subsequently we extend the con- 
cept of optical SUSY to two-dimensional (2D) settings 
with cylindrical symmetry, as in optical fibers. We show 
that a partner potential with a SUSY spectrum of radial 
modes exists, offering the possibility for angular momen- 
tum multiplexing. Finally, we investigate the implica- 
tions of supersymmetry within the framework of finite 
periodic structures and propose a versatile approach to 
systematically design SUSY optical lattices. 

To explore the consequences of supersymmetry in op- 
tics, we consider optical wave propagation in an arbi- 
trary one-dimensional refractive index distribution n{x). 
Waves propagating in the xz-plane can always be de- 
composed in their transverse electric (TE) and trans- 
verse magnetic (TM) components. For TE waves the 
field evolution is governed by the Helmholtz equation 
{dxx + dzz + kQn'^{x)) Ey{x^ z) = 0. Modes propagating 
in this system have the form Ey{x^ z) = f{x)e^^^ and sat- 
isfy the following eigenvalue equation for the propagation 



constant /3: 



nf{x) = -p^fix), 



(1) 



where H = — ^ — kQin?{x) corresponds to the Hamil- 
tonian operator in a Schrodinger equation. For a given 
index profile rS^\x)^ SUSY now provides a systematic 
way for generating a superpartner n^'^\x). If the index 
distribution n^^\x) supports at least one bound state 
/} ^ (x) (the ground state) with a propagation eigenvalue 

/3f \ SUSY can be established via H^^^ + (p^^^\ = At A, 

where A = -^d/dx + W{x) and A^ = —d/dx + W{x) are 
defined in terms of a yet to be determined superpotential 
W{x). The optical potential and its superpartner then 
satisfy 



(W''H^))' = (/^i'^)' - ^' ± ^'. 



(2) 



Taking into account that A^Af^ =0, one finds that 
^/i — 0. Thus a valid solution for W can be obtained 
from the logarithmic derivative of the node-free funda- 
mental mode: 



W{x) 



dx 



In 



(fi'Hx)). 



(3) 



Figure [TJa) depicts an arbitrary refractive index dis- 
tribution supporting a set of six guided modes. Here 
the maximum index contrast is 5 x 10~^ and the wave- 
length used is 1pm. While Eqs. (1-3) are valid in the 
Helmholtz regime, here we consider a low contrast struc- 
ture that is experimentally feasible. For this example, 
the SUSY partner (Fig. [ijb)) has been numerically cal- 
culated from Eq. ([2| through the corresponding super- 
potential (Fig. Illc)) that was obtained from Eq. ([3|. As 
Fig. ma) clearly shows, the fundamental mode of n^^^ 
lacks a partner in the eigenvalue spectrum of n^^\ indi- 
cating unbroken SUSY. On the other hand the second 
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FIG. 1. (Color online), (a) Exemplary refractive index land- 
scape (grey area) and its six bound modes (vertical placement 
indicates their respective eigenvalues), (b) SUSY partner and 
its five modes. The operators A, A^ transform the phase- 
matched modes into each other, (c) Both index landscapes 
can be constructed from the superpotential W and its slope 




Multimode Lossy 
Waveguide Superpartner 



FIG. 2. (Color online). Beam propagation in a multimode 
waveguide, (a) When isolated (before dashed line), and when 
coupled to its lossy superpartner (after dashed line, losses: 
a ~ 0.4cm~^). Two more advanced stages of this same field 
evolution in the coupled system are shown in (b,c). 



state of n^^^ is paired with the first mode of n^^^ that 
has exactly the same propagation constant in spite of 
its different parity. In this way, all the modes of these 
two superpartners can be perfectly phase-matched ex- 
cept for the fundamental mode of n*^^^. Therefore SUSY 
provides the only strategy we know of to achieve global 
phase matching conditions, irrespective of how large the 
number of modes is, in such multimode optical poten- 
tials. 

This latter feature can be exploited for mode filtering 
applications. The idea is illustrated in Fig. |2Va) where 
n^^^ has the form of a step-index like waveguide that 
supports three modes at A = 1.5pm. The optical propa- 
gation when this system is excited with an arbitrary in- 
put beam, is depicted in the first propagation section of 
this figure. In this range, the field evolution is seemingly 
chaotic because of modal interference. Once however the 
superpartner waveguide is put in proximity, all the modes 
of n*^^^ (apart from the fundamental) are periodically cou- 
pled between these two structures. Despite their par- 
ity, coupling between the phase-matched modes occurs 
through their evanescent tails. If for example the sec- 
ond waveguide is made intentionally lossy, all the modes 
of n^^^ eventually disappear except the fundamental, as 
shown in Figs. |2|b,c). Similarly, the fundamental mode 
can be selectively amplified. This behavior could be po- 
tentially useful in large mode area laser sources. 

SUSY structures also exhibit identical scattering prop- 
erties in terms of their refiectivities and transmittivities. 
In this case, the radiation mode continua are related to 
each other through the SUSY algebra. Let us consider 
again the SUSY pair described by Eqs. ([2|. We also 
assume that n^^^ asymptotically approaches a constant 
background value Uoo at x -^ ±00. For an angle of inci- 
dence ^, the components of the incident wave vector are 
kx = koTioo cos{6) and kz = A:o^oo sin(^). The SUSY for- 
malism then relates the field reflection/transmission co- 
efficients r^^'^^ and t^^'^^ associated with these two struc- 
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FIG. 3. (Color online). Scattering properties of the SUSY 
pair from Fig. IT] (a) logarithmic plot of the angle- dependent 
reflectivities EP^''^-' (the graphs have been offset for visibility), 
and (b) Phase difference of the transmission coefficients t^^'^-* 
(inset: Schematic of the scattering configuration) 



tures in the following way [15]: 



.(2) 



where Wr^ 



W^ 



±i^^(i) ^(2) 



Woo + iK 



w^ 



ikn 



t<«, (4) 



(tf')- 



2^2 



k^n 



represents the limit of 

the superpotential at x ^ +oo as obtained from Eqs. ([2|. 
Note that the argument of the square root is always a 
non-negative quantity [16]. It follows that n^^'^^ exhibit 



i?(2) 



^^-'■'^^ I and trans- 



identical reflectivities i?^^^ 

mittivities T^^) = T^^) = [t^^'^^l^ Consequently, bar- 
ring direct phase measurements, the two SUSY structures 
would be indistinguishable at any angle of incidence. In- 
terestingly, the phase difference between r^^^ and r^^\ 
and between t^^^ and t^^^ for any given is solely deter- 
mined by the propagation constant /3} ^ of the fundamen- 
tal mode and the background refractive index Uoq. 
A schematic of a possible scattering arrangement is 
depicted in Fig. |3] The angle-dependent reflec- 
tion/transmission coefficients for the SUSY pair consid- 
ered in Fig. [l](a,b) were evaluated by means of the differ- 
ential transfer matrix method [17 when the background 
refractive index is rioo = 1.5. In accordance with our 
previous discussion, the two structures display identical 
reflectivities (Fig. [sj^a)). The phase difference between 
their respective transmission coefficients is also shown in 
Fig-ib). 

Having investigated SUSY in ID optical systems, the 
question naturally arises as to whether these concepts can 
be extended to 2D structures. The answer is not particu- 
larly obvious given that the aforementioned factorization 
technique relies on ID Hamiltonians [8 . In what fol- 
lows, we show that this limitation can be overcome in 
paraxial settings with cylindrical symmetry, as in weakly 
guiding optical fibers. In this regard, let us consider the 
radial refractive index profile n{r) = Uoo + An(r) where 
An <C rioQ. In this case, the slowly varying field envelope 
U satisfies the paraxial equation 



_9^_ 1_9__^_9^_ 

^^2 ^ Qrjj ^2 g^2 



u 



4^' 



(5) 



where r] = r/ro is a normalized radial coordinate, r^ is an 
arbitrary spatial scale, (j) is the azimuthal angle and the 
normalized axial coordinate is given by ^ = z/(2/conoo^o)- 
In this representation the optical potential reads V = 
^nook^rQAn. By expressing the mode U = e^^^e^^^R{r]) 
in terms of its orbital angular momentum ^, and after 
using the radial transformation R 
Eq. ([5| to a ID form. 



T] ^I'^u we reduce 



-^-^eff(r?))« 



-\xu, 



(6) 



with the effective potential Kff(7^) = y{r\) H — —^ — • 
By designating the modes of Eq. (|6| as ugm^ having 
azimuthal and radial mode numbers i and m respec- 
tively, one can then generate an effective partner poten- 
tial V^^ (r\) for a given effective potential V^^ (77). As 
in the ID case investigated before, these two potentials 



,(1) 



are related via the fundamental mode u\ \ of the first 

potential; V^ = V^ + 2^ In f'^^^i). In the original 

coordinate system, Iv^ \ = r]~^^'^u\ {, which yields the 
following relation between the superpartner potentials: 



V^^\r]) = V^^\r]) 






-R 



(7) 



Note that in deriving the most general expression for F^^^ 
we have assumed a different azimuthal mode number £2 
for the partner potential. In other words, a potential 
V^^^ and its partner V^^\ constructed for a certain ii 
and ^2, will only be supersymmetric with respect to the 



?(i) 



?(2) 



subsets R\ ^ and R\ ^ of their respective radial modes 
(m = 1,2,...). Note that the second term in Eq. ^ 
may introduce a singularity at 77 = 0. Yet, this can be 
alleviated through an appropriate choice of ii and £2- 
Near the origin {rj <C 1), the radial solutions Rirn of 
any well-behaved potential V{r]) are proportional to 77 1 ^1 
[15], and thus Ri,i{r]) ^ ^'^''- Therefore, Eq. ^ yields 
a non-singular partner potential only if |£2| = Ki| + 1- 
This relation reveals an unexpected result; in cylindri- 
cally symmetric settings, SUSY provides a link between 
sets of modes with adjacent azimuthal numbers. Given 
that V^^^ vanishes at 77 ^ 00 it then follows that [15] 
^exp {-^JJTlari), and hence V^'^\r]) ^ l/rf in 



(1) 



this same limit. 

Figures |4]^a,b) depict the field profiles of the modes 
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(1,2) 



g*^i,20^^ ' ^(77) corresponding to the two 
cylindrical superpartner index profiles in Figs. |4|^c,d). 
In this case, the original refractive index distribution is 
taken to be /S.n(r) = (5e~^^/^°^ , where the core radius is 
ro = 30pm, the index contrast amounts to S = 2 x 10~^ 
and the background refractive index is Uoo = 1.5. At 
a wavelength of 1.55pm, it supports a total of twelve 
guided modes. Based on the lowest state with ^1 = 1, a 
partner potential for ^2=2 was generated according to 
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FIG. 4. (Color online). (a,b) Supersymmetric subsets of 
bound states corresponding to the SUSY pair of cylindrically 
symmetric index profiles (c,d) generated for azimuthal num- 
bers ^1 = 1 and ^2 = 2. (e,f) Complete eigenvalue spectra 
(effective refractive indices) of both potentials. The respec- 
tive subsets of SUSY states are indicated by dashed frames. 



Eq. [see Fig.gd)]. Note that whereas SUSY holds 
between the modes with ii = I and -^2 = 2, the rest of 
the eigenvalues remain disjoint, as shown in Figs. [4Fe,f). 
By relating mode subsets of different azimuthal indices 
in this 2D setting, SUSY offers the possibility for a fully 
integrated realization of optical angular momentum mul- 
tiplexing |18j. 

We next consider SUSY in finite periodic arrangements. 
For example, a lattice of N well-separated single- mode 
waveguides is known to support a set of N bound states 
or supermodes. In this array environment, the funda- 
mental state is again node-free and hence can be readily 
used to generate a superpotential according to Eqs. ^. 
The corresponding SUSY partner resembles a lattice with 
N — 1 dissimilar channels located in the gaps between the 
original waveguides (see [15 ). 

The coupled mode formalism provides an effective way 
to describe wave evolution in photonic lattices within the 
first band. The set of coupled differential equations [19] 
for the modal field amplitudes a can be written in the 
form Hs. = Aa, where H is now the discrete Hamiltonian 
of the system. This discretization provides a powerful 
approach for constructing SUSY pairs: The Hamiltonian 
can be directly factorized using the Cholesky method 
[20] . The pair of isospectral Hamiltonians thus obtained 
retains the tri-diagonal shape of H, i.e. the SUSY partner 
represents again a photonic lattice with nearest-neighbor 



coupling. Note that whereas both Hamiltonians are 
N X N matrices, SUSY is nevertheless unbroken in the 
sense that the N^^ waveguide of lattice 2 is completely 
decoupled. 

Even more importantly, the discrete formalism outlined 
above relaxes the need for exactly controlling the refrac- 
tive index landscape. In particular, the technological 
difficulties associated with sharp index depressions can 
be circumvented without any loss of functionality. In- 
deed, the control of only two parameters is here sufficient 
for the actual realization of SUSY optical systems: The 
waveguide's effective refractive index, which determines 
the propagation constant, and their separation, which 
relates to the coupling coefficient. A sequence of SUSY 
potentials can be iteratively obtained by discarding the 
respective isolated channels. Such a SUSY "ladder" can 
facilitate a lossless decomposition of any input beam into 
its modal constituents. A weak coupling cl between such 
consecutive partner lattices, as indicated in Fig. [sja), 
does not perturb SUSY and allows for an interaction only 
between states with equal eigenvalues. Consequently, en- 
ergy initially carried by the k^^ supermode in the funda- 
mental lattice can be transported between all layers l.../c, 
but is rejected by layer k-\-l. The propagation dynamics 
arising from the excitation of several supermodes in the 
fundamental lattice are shown in Figs, [sfb-d) for such 
a SUSY ladder based on a uniform array with Nq = 6 
waveguides. The condition of weak inter-layer coupling 
was assured by setting cl to be 5% of the coupling C 
within the uniform lattice. 

In conclusion we have shown that SUSY partner systems 
can be generated for any ID refractive index landscape 
supporting at least one bound state. Despite their dis- 
similar shapes, SUSY structures can exhibit identical re- 
flectivities and transmittivities for arbitrary angles of in- 
cidence. Subsequently the concept of optical SUSY was 
extended to 2D settings with cylindrical symmetry. In 
this case SUSY was established for sets of modes ex- 
hibiting consecutive azimuthal indices. In the context 
of photonic lattices, SUSY manifests itself as a reduction 
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FIG. 5. (Color online), (a) Schematic of a SUSY ladder with 
N = 6 layers. Propagation dynamics when a supermode of 
the original lattice is selectively excited, (b) /c = 1 (fundamen- 
tal state): Confined in the first layer; (c) A: = 3: Penetrates 
only the first 3 layers (d) /c = 6: Moves freely across the entire 
ladder 



in the number of channels. This concept is general and 
highlights the potential of SUSY for robust optical filter- 
ing and signal processing applications. 
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